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Abstract 

Amongst the class of supergravity solutions found by Lin, Lunin and Maldacena, 
we consider pure and mixed state configurations generated by phase space densi- 
ties in the dual fermionic picture. A one-to-one map is constructed between the 
phase space densities and piecewise monotonic curves, which generalize the Young 
diagrams corresponding to pure states. Within the fermionic phase space picture, a 
microscopic formula for the entropy of mixed states is proposed. Considering ther- 
mal ensembles, agreement is found between the thermodynamic and the proposed 
microscopic entropies. Furthermore, we study fluctuations in thermodynamic en- 
sembles for the superstar and compare the entropy of these ensembles with the area 
of stretched horizons predicted by the mean fluctuation size. 
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1 Introduction 



Extremal small black holes are intriguing objects. On the one hand, as solutions of super- 
gravity (SUGRA), they have vanishing horizon area, because the horizon locus coincides 
with the space-time singularity. On the other hand, counting microstates yields a finite 
entropy (for a review, see [1]). Within String Theory, the origin of this discrepancy is to 
be found in the supergravity approximation. In simple terms, in the supergravity limit, 
one considers objects of the size comparable to a typical gravitational length scale, L, 
which is much larger than the string length, L ^> a' 1 ^ 2 . Consequently, objects of a typical 
size A<X cannot appear, even if A ^> a . This is the case of small black holes. The 
horizon area that would follow from their entropy and the Bekenstein-Hawking formula, 
A = AGS, vanishes in the supergravity limit when measured in units of L. Including 
higher derivative corrections to SUGRA generates intermediate length scales; e.g., one 
could have A 4 = a'L 2 . Therefore, for small black holes, a finite-size horizon is expected 
to be generated, and the coefficient 1/4 in the Bekenstein-Hawking formula can receive 
large corrections [2]. 

Maybe the easiest way conceptually to find the entropy of extremal small black holes is 
to use Sen's entropy function [3] for a SUGRA Lagrangian with higher derivative correc- 
tions pfl[5] El Another way, which yields, however, only an order-of-magnitude estimate of 
the entropy, employs the fuzzball conjecture (for a review, see [9]). According to this con- 
jecture, the black hole geometry is the result of coarse graining the underlying microstates, 
which possess distinct geometries differing from each other only in the space-time region 
around the singularity. To observers far from the black hole measuring typical observables 
(in the sense discussed in [ID]), the differences between microstates are invisible. Placing a 
stretched horizon around the space-time region, where the microstates significantly differ, 
one can estimate the black hole entropy using the Bekenstein-Hawking area law. 

In this article, we consider the bubbling solutions found by Lin, Lunin and Maldacena 
(LLM) which are the gravity duals of 1/2-BPS chiral primary operators in M = A 
Super Yang-Mills (SYM) theory. The 1/2-BPS sector has a simple description in terms of 
free fermions in a harmonic oscillator potential, and an orthogonal set of states is in one- 
to-one correspondence with Young diagrams [T2"j [T5] . The LLM bubbling solutions have 
been extensively studied in the literature. Duality between the Pauli exclusion principle 
and the absence of closed time-like curves has been established in [HI [T5] . Statistical 
ensembles of bubbling solutions have been considered in [T6l [TFl [TEl [TO]. Further work 
can be found in [Hll2niI2IlI22lI231I21lI2SlI2niI2Z]- 

In [TU] , an ensemble of microstates has been proposed that fits the geometry of the su- 
perstar [2HI [29l [30] . In the present article, we will further test this ensemble by calculating 
an entropy estimate using a stretched horizon, lifting the geometric analysis of [32] to ten 
dimensions. To this end, we shall consider the mean size of fluctuations of the superstar 

X A similar method is c- function extremization [B]. In this method, higher derivative corrections have 
been considered in [7]- A generalization of Sen's entropy function formalism for rotating black holes has 
been developed in [5]. 
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ensemble, which provides a measure for the size of the stretched horizon, with the result 
that the predicted horizon size is too large in orders of magnitude. Hence, we will propose 
another ensemble, which we shall call the restricted superstar ensemble, which also gives 
rise to the superstar geometry, but whose fluctuation sizes are consistent with the order 
of magnitude of the superstar entropy. On the way, we will discuss in detail the map be- 
tween Young diagrams that characterize microstates and fermion droplets that generate 
the LLM geometry, as well as the generalization of that map to mixed states. Moreover, 
we will propose and test a unique formula for the entropy of mixed states, which arises 
directly from identifying the LLM plane of droplets with the semi-classical phase space 
of fermions. 

The article is organized as follows. In Sec. we briefly review the LLM solutions and 
some of their properties. Sec. [H] deals with the one-to-one map between fermion phase 
space density functions u(r) and the (piecewise) monotonic curves Y(X), which generalize 
the Young diagrams associated with microstates. At the end of that section, we propose 
our new entropy formula. In Sec. HI we discuss the construction of a stretched horizon in 
the superstar geometry, lifing the 5d analysis of [TTJ to ten dimensions. This will allow us 
to make contact later with the typical size of fluctuations in ensembles of Young diagrams. 
In Sec. El some statistical ensembles of Young diagrams are considered. We shall start 
with a review of the general grand canonical ensemble of [10J, by means of which we will 
be able to test our entropy formula. Then, we shall specialize to the grand canonical 
superstar ensemble of [10] and find that the size of the stretched horizon compatible with 
the fluctuations of that ensemble is too large. Therefore, we define and study a restricted 
superstar ensemble, whose typical fluctuation size matches the expected size of a stretched 
horizon. Finally, Sec. [6] will contain conclusions. 

2 Review of LLM solution 

The metric of the LLM solutions is given by 

ds 2 = -h- 2 {dt + Vidx*) 2 + h 2 {dy 2 + dxW) + ye G dQ 2 3 + ye~ G dQ 2 m 
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with i — 1, 2, and 



h- 2 = , 2y , e 2G = i±^. (2) 
VI - 4z 2 ' l-2z 1 ; 



The auxiliary function z and the vector Vi are given by 

zfcv) = — d x 



n J [(x-x') 2 + y 2 ] 2 ' 

7T J [(x — x')^ + y z \ z 

where zq(x) = z(x, 0) is the boundary data for the auxiliary function z(x, y) on the plane 
y — 0, spanned by x 1 and x 2 . For the complete solution including the 5- form, we refer 
the reader to the original paper [TT] . 
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The metric (CQ) is regular, if zo( x ) = ±1/2. It possesses naked null-singularities for 
\z 1 < 1/2, and closed time-like curves (CTCs) for \z \ > 1/2 [14j@ In this paper, we shall 
consider only solutions with \z \ < 1/2. 

It is convenient to define ^ 

zo(x) = - - u(x) , (4) 

so that, for a regular metric, only w(x) = 0, 1 are allowed values. For better visualization, 
one can think of u as specifying a certain darkness, or grayscale, on the 1-2-plane, with 
u = and u = 1 corresponding to colours white and black, respectively. Then, the regular 
solutions are characterized by "droplets", i.e., arbitrary configurations of black and white 
regions on the 1-2-plane, whereas gray areas generate naked null-singularities. 

Intuitively, one can identify it(x) with the semi-classical phase-space density of fermions 
in the dual fermionic picture that describes the 1/2-BPS sector. In this context, the Pauli 
exclusion principle for fermions (and holes), < u < 1, implies the absence of CTCs |14j . 
The regular solutions described by a distribution of black droplets are the duals of pure 
states, or microstates, whereas solutions that involve gray areas are duals of mixed states. 

In what follows, we shall consider solutions corresponding to a large, but finite, number 
of fermions, N. Identifing the 1-2-plane with the semi-classical phase space of the fermions, 
we should remember that it is "discretized" into small elements of volume 2nh. For small 
h, we approximate the sum over these phase space elements as 

1 



E 



2nh 

ph.sp. 



d 2 x . (5) 



Therefore, the total number of fermions is 

1 

2nh 



E 



u = — -r / d^rw(x) = N . (6) 



ph.sp. 

Without loss of generality, the origin in the 1-2-plane is identified with the center of mass 
of the density u(x). 

For large y 2 + |x| 2 , the metric (pQ) approaches asymptotically AdS 5 x S 5 , with length 
scale 

'4 _ 1 / j2 



L=^j d J xu(x) . (7) 

Hence, a circular black droplet of radius R corresponding to the fermionic ground state 
(the Fermi sea) generates precisely AdSs x S 5 , with L 2 = R. 
From dHD and (J7j) one readily obtains 

L 4 = 2HN . (8) 

Comparing this with the standard formula, L 4 = 4ng s a' 2 N, yields the "Planck constant" 

k = 27rg s a' 2 . (9) 



2 Clearly, one must cut out those regions of the (x, y) 3-space, where \z(x, y)\ > 1/2. Then, the 
metric is singular at the space-time boundary, where \z(x, y)\ = 1/2 for y > 0. CTCs are found within 
space-time, close to this boundary. 
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The leading deviations from the asymptotic AdSs x S 5 geometry reveal the presence of 
equal angular momentum and mass [11], given by 

A = 3 = ss / d2 ^ M(x) - Ks« / ^ u(x) ) 2 ■ (10) 

This corresponds precisely to the energy of the fermions above the ground state. 

Finally, considering regular solutions, the quantization of 5-form flux imposes an ad- 
ditional restriction on the droplet distribution. It can be shown that the area of any 
compact black or white region is quantized in Planck units [TT] . 



area) = 2-nhn , (11) 



with some integer n. 



3 Droplets, Young diagrams and limit curves 

In this section, we shall consider LLM solutions determined by density functions u(r) with 
rotational symmetry on the 1-2-plane. Such solutions possess two Killing vectors and are 
static in the asymptotic AdS coordinates. We will show that there exists a one-to-one 
map between these LLM solutions and (piecewise) monotonic curves X(Y). For regular 
geometries generated by distributions of concentric black rings, these curves precisely 
delimit the Young diagrams that specify the dual 1/2-BPS microstatesl For general 
densities u(r), they are the limit curves of [TU] . 

In what follows, we use polar coordinates (r, 0) on the 1-2-plane. 

3.1 Microstates: Droplets and Young diagrams 

Let us start with a droplet configuration of M concentric black rings, such as those 
illustrated in Fig. [TJ Since the LLM geometry generated by such a configuration is regular, 
it is interpreted as a microstate, i.e., a pure quantum state. The quantization of five-form 
flux (TTTT) introduces two sets of integers, b k and w k (k = 1, . . . , M), specifying the areas 
of the single black and white rings, respectively. More precisely, if we denote by r k and 
Rk the inner and outer radii of the fc-th black ring, respectively, where k — 1 stands for 
the innermost ring, then the integers b k and w k are defined by 

h = ^(Rt ~ rl) , w k = JL*(r2 - RU) . (12) 

Here and henceforth, we define Rq = 0. Moreover, r\ = is an allowed value corresponding 
to a configuration where the innermost ring is a disc, as in Fig. [T^,). The total area of all 
black rings specifies the total number of fermions, so 

M 

Y^b k = N. (13) 
fe=i 

3 In this case, the map can be found explicitly in Fig. 1 of 19J. 
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a) 



X 




Figure V. Examples of black ring distributions and their associated Young diagram curves. 

Using the integers bk and Wk, one can construct a curve that delimits a Young diagram 
as follows. Take a pen and place it at the origin of a usual X-Y coordinate system^ 
Without lifting the pen, draw a line Wi units to the right, followed by 61 units to the top, 
then W2 units to the right and 62 units to the top, and so on for all k = 1, ... , M. From 
where the pen is then, draw a horizontal line extending infinitely to the right. This line, 
which is at Y = N (in as yet unspecified units), corresponds to the infinite white plane 
around the droplets. The area between the curve, the V-axis and the horizontal line at 
Y = N (gray areas in Fig. [T]) constitutes the Young diagram that specifies the 1/2-BPS 
operator dual to the LLM geometry. Fig. [1] shows two examples of this construction. 

In terms of the integers bk and Wk, the area of the Young diagram is 

M k 

A = 5>5><- ( 14 ) 

k=l 1=1 

4 We use capital letters X and Y in order to avoid confusion with the coordinates in the LLM solution. 
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Let us compare this formula with LLM's formula ( jTUi) . From ( jTUl) we obtain 



k=l 
M 



M 



2 



k=l 



k=l 
M 



(Rl + rl)-(Rl-rl)-2j2(Rt-r?) 



Kk 



4ft 2 



(15) 



which coincides with (1141) after using (112p . 



3.2 General states: Fermion densities and limit curves 

It is possible to generalize the above construction to the case of rotationally symmetric 
phase-space densities u(r). For this purpose, we follow the proposal of [TU] and define two 
functions Y(r) and X(r) bw\ 



Y(r) = J dr'r'u{r') , X{r) = K 2 - Y{r) . (16) 
o 

Notice that both Y(r) and X(r) are monotonically increasing, because < u < 1. The 
maximum value of Y is 

Yoo = Y(r = oo) = hN = -L 4 . (17) 

The pair [X(r), Y(r)] defines a parametric curve on the X-Y plane. It is easily seen that, 
for a distribution of concentric black rings such as we considered earlier, ( ITBT) reduces to 
the construction of the Young diagram curve considered in the last subsection, with unit 
length h on both axes. Fig. [2] shows a simple example containing areas with u — 1/2. 

We can also go the other way around and find the fermion density uir) associated 
with a given monotonic curve Y(X). To do this, we invert ffTB"]) assuming that X is a 
(piecewise different iable) function of Y. Then, one obtains 

r=y/2(X + Y), u = 1+ 1 Xf(Y) • ( 18 ) 

Similarly, one can write 

t—T+hm- (19) 

where Y is considered as a (piecewise differentiable) function of X. 

Expressing the LLM energy (imp in terms of Y and X using (1161) and (1181) . one obtains 

A = ±JdYX(Y), (20) 



5 With respect to [10j . we interchange X and Y and use a different normalization. 
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Y 




X 

Figure 2: The grayscale distribution obtained by superposing the distributions of Fig. [T] 
with equal weights, and its associated curve. 



which is, as expected, the area to the left of the curve Y(X) measured in units of h. 

Now, we would like to sharpen the notion that singular geometries describe mixtures 
of microstates. To this end, we first establish that the above construction satisfies a nice 
superposition property. In the following, let h and iV (and, consequently, also L) be fixed. 
Consider a number of rotationally symmetric density functions u k (r), each describing a 
phase space distribution of N fermions. These u k can correspond to microstates, but they 
need not. Let us now consider the following fermion density, 



a k u k (r) 



with a k = 1 



(21) 



Then, from (1161) follow immediately 



Y(r) = J2a k Y k ( 



(22) 



i.e., the superposition of fermion densities is equivalent to the superposition of Young 
diagram limit curves, if the latter are regarded as parametric curves depending on r. 

It is worth noting that this superposition property is consistent with the notion of 
microstates as purely black and white droplet configurations. Indeed, ( 1211) implies that it 
is impossible to obtain such configurations as superpositions. 

In order to show that the energy A associated with the mixed state with density u(r) 
is given by the weighted average of the energies A^, one can start from (I20I) and use ([161) 



to rewrite it at 



A 



h 2 



dY -r — Y 




drr u(r) — 



2^Z) Q fe f J drr 3 u k {r)-Yl \ = ^ a k A k . 



(23) 



6 Notice that X(Y) ^ J2k a kX k (Y). 
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3.3 Entropy formula 



Mixed states possess a non-zero entropy. Therefore, it should be possible to find a unique 
measure for the entropy of a mixed state, given its density function u(r). The von Neu- 
mann entropy of mixing, S m ix = — J2k a k l na fc ; where one assumes that the densities Uk in 
([21]) correspond to pure states, is not a suitable definition, however. One reason is that, 
given just the mixed state density u(r), it is impossible to determine unambiguously its 
microstate content and the coefficients a^. 

A good definition for the mixed state entropy follows from the interpretation of the 
droplet plane as the phase space of fermions. Remember that, semi-classically, the phase 
space is quantized in small elements of volume 2nh. Let us label these elements by an 
integer n. Then, the probabilities for finding or not finding a fermion in the phase space 
element j^n form a Bernoulli distribution, {u n , 1 — u n }. Summing the von Neumann 
entropies of the Bernoulli distributions of all phase space elements, we obtain an unam- 
biguous entropy for the mixed state configuration. In the semi-classical limit, this gives 

S = -—r [ d 2 x[u\nu + (1 - «) ln(l - «)] . (24) 
27m J 

We shall confirm in Sec. [5] that f [24"l) correctly reproduces the thermodynamic entropy for 
ensembles of Young diagrams. 



4 Superstar and stretched horizon 

We will consider now the construction of a stretched horizon for the superstar geometry, 
which provides an estimate for the entropy of the underlying mixed state. Let us start 
with the 5d reduction of the superstar, which is a well-known extremal black hole solution 
of M = 2 gauged SUGRA in five dimensions [28, 29j. Its metric is given by 

dsj = -H^fdt 2 + H 1 / 3 (r Mp 2 + p 2 d{lj) , (25) 

where 

T 2 2 

H(p) = l + u J —, f{p) = l + H{p)^. (26) 

For convenience, we have expressed the black hole charge by the dimensionless u in units 
of the AdS length scale, L. We use p instead of the conventional r in order to avoid 
confusion with the polar coordinate r on the 1-2-plane in the lOd solution. 

Following [T7j, let us consider a static stretched horizon located at constant p = p 
and wrapping the 3-sphere in ( 1251) . Its "area" is easily found to be 

A 3 = 2n 2 p 3 H(p y/ 2 . (27) 

Because the 5d J\f = 2 gauged SUGRA can be thought of as the compactification of lOd 
type-IIB SUGRA on a 5-sphere with radius L, the 5d gravitational constant is 

G W _ 2^h 2 _ vr 
G "-UV^- LV -2W L ' (28) 
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where we have made use of (jHD and (jHJ). Hence, an estimate for the entropy is found as 



A. = 7r jv 2 ^y^+^. (29) 

It is instructive to repeat this calculation in ten dimensions using the LLM form of 
the superstar solution. After uplifting the 5d metric ( l25l) . one obtains the lOd superstar 
in the form 130 



ds l = __L ( cos 2 e + D ^) dt 2 + lL sin 2 edtd(j) + EE r sin 2 ^ d0 2 + 

+ v / ^D(/- 1 dp 2 + p 2 d^) +L 2 v / ^d^ + 4^cos 2 0d^ , 



where D = sin 2 9 + H cos 2 9, 9 G [0, 7r/2], and G [0,27r]. Upon changing coordinates to 

M 

y = Lp cos 9 , r = L 2 a/ f(p) sin 9 , t — > , (31) 

one finds the superstar in the standard LLM form (JTJ). The 1-2-plane is described in polar 
coordinates (r,<p). The fermion distribution generating the solution is given by a gray 
disc of radius 

r = L 2 y/l+u (32) 

and uniform density 

u = {l+u)~ 1 . (33) 

Hence, according to fllSp . the limit curve corresponding to the superstar is a straight line 
with slope X'(Y) = u, for < Y < Y^, and the delimited area is just a right angled 
triangle with side lengths = hN and X(Foo) = hNui. The energy of the superstar 
follows easily from ( 1201) . 

A = ^N 2 u . (34) 

According to our general formula (EMI) , the entropy of the underlying mixed state is 

S = JV[(1+ w) ln(l + w) -u\nu] . (35) 

Let us consider now the uplift of the stretched horizon. It is a static 8-surface wrapping 
the two 3-spheres and forming a rotationally symmetric 2-surface above the y — plane. 
Its y-r profile follows from (IBTl) and is the segment of an ellipse, 

v 2 r 2 

\ + — = l } (36) 

VH r H 

with 

y H = Lp , r H = L 2 ^Jl + u + p 2 /^ . (37) 
From the LLM metric ([T]) one obtains the area of the stretched horizon as 



As = J an 3 an 3 d</> dr y 2 y/i + ^J f~ - z 2 ) r 2 - y w 2 , (3? 
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where y' = dy/dr. To continue, we still need expressions for z and on the horizon. 
The easiest way to find them is from the uplift formulae [30] , which imply 

1 p 2 D - L 2 cos 2 6 l pt + (u- l)y 2 
Z 2 p 2 D + L 2 cos 2 9 2 pi + (u + l)y 2 ' 1 J 

cos 2 + p 2 D/L p A + {uJ + l)y 2 ' 1 1 

where the second equalities in each line follow from (15T1) . After substituting these formulae 
into (1551) . changing integration variable from r to y and using ( 1561) and ( 1571) . the expression 
simplifies dramatically after some algebra. We end up with 



VH 2 



As = (27r 2 ) 2 (2n) [ dyy^Ju + p 2 /L 2 = 2^ p 2 Q L & J uj + p 2 /L 2 



Po 



Hence, we obtain for the entropy 



(41) 



27r 5 p 2 L + p 2 /L 2 r2 p 2 
hor ~ 4G 10 ~ 4(2vr^2) - ^ L 2 



yf^+pJ/L 2 , (42) 



in agreement with ( 1291) . 

As we have seen in ( 1351) . the entropy is of order N. In order to achieve this, the radius 
parameter p must be of order 1/y/N in units of the AdS length scale L. More, precisely, 
setting 

<«' 

where c(u) can depend on the slope of the limit curve, and neglecting the second term in 
the square root (we have large N), we obtain 

S hor w c(uj)N^ = c(cj)v / 2A , (44) 



where A is the superstar energy (J5J 

One can justify the form (1431) by considering the minimum fluctuation of the number 
of columns of a Young diagram [T7j, which is, obviously, 5n c = 1. Fluctuations of the 
number of columns translate into fluctuations at the upper end of the limit curve triangle. 
Moreover, as the upper end of the limit curve corresponds to the border of the gray disc of 
non-zero fermion density, we can interpret these fluctuations as fluctuations of the disc's 
radius. In particular, let us use (ITU]) to writdj] 

5r 2 = 2 5X(Y O0 ) = 2h 5n c . (45) 

Considering a stretched horizon of the form (15U|) . we identify 

r% = rl + Sr 2 . (46) 

7 Y nc , is fixed. 
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Thus, using also (152"]) and ( 1371) . we determine the parameter p Q as 

§ = T^ r «-^ = T^-Jf Sn " (47) 

implying c(u>) = 7r in (HHI) and (IHJ). 

Using the minimum fluctuation of a statistical ensemble is, however, not very rea- 
sonable. In Sec. El we will check whether the mean size of fluctuations in the superstar 
ensemble is in agremment with (H31) . 



5 Ensembles of Young diagrams 

In this section, we shall consider some statistical (thermodynamic) ensembles of Young 
diagrams. Our main interest is in ensembles, whose limit curves are in agreement with the 
superstar solution. We will, however, start by considering the grand canonical ensemble 
of [TO] and verify that the fermionic entropy f[2~4"j) reproduces the thermodynamic entropy 
of the ensemble. 

First, let us state how we describe Young diagrams. The most convenient way for our 
purposes to describe a given Young diagram is to specify N numbers Cj, j = 1, 2, . . . , N, 
which count the number of columns of length j. The maximum possible length of any 
column is N. The total number of columns of the diagram thus specified is 

N 

n c = Y, c i> ( 48 ) 
i=i 

whereas the total number of boxes of the diagram is 

N 

E = J2jCj. (49) 
i=i 

For 1/2-BPS states in M = 4 SYM with gauge group U(N), the c/s can take any non- 
negative integer value, Cj = 0, 1, 2, . . .. For gauge group SU(N), the restriction tr Z = C§ 
removes all diagrams with c\ > 0. In what follows, we shall consider U(N). In the 
large- N (thermodynamic) limit, the differences between the U(N) and SU(N) results are 
sub leading in 1/N. 

In the thermodynamic limit, we replace the integers j and c, by the continuous vari- 
ables of the limit curve, such that 

j 

hj^Y, hY.^^XiY) Cj ^X'(Y), (50) 

i=i 

and take the limit such that 

h -> , iV -> oo , with hN = fixed. (51) 
The last replacement in ( 1501) comes very handy in view of the relation (fl8l) . 

S Z is the complex combination of two of the six real scalar fields, e.g.,Z — X 1 + iX 4 , which transform 
in the adjoint of the gauge group. 
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5.1 Grand canonical ensemble 

Balasubramanian et al. [ID) studied a grand canonical ensemble, in which the mean energy, 
A = (E), and the mean number of colums, N c = (n c ), were fixed using Lagrange multi- 
pliers. They found that the superstar limit curve is obtained for infinite "temperature" . 
In the following, we shall briefly review their calculations for arbitrary temperature, in 
order to verify our proposal for the entropy fl24"l) . and specialize to infinite temperature in 
Sec. El 

Let us start with the grand canonical partition functional 

N 

■^-^ , ^-^ - E (aj/N+\)ci 

Z = Y, e- aE ' N ~ Xn - = e h ■ (52) 

{ Cj } { Cj } 

After performing the sum over the c/s in (|52|) . we obtai 

N 

Z = X{[l-e- {a3lN+X) Y 1 ■ (53) 

To continue, we consider InZ and take the thermodynamic limit ( 15"T|) . with the result 

lnZ = - [Li 2 (e~ A ) - Li 2 (e-( Q+A ))] , (54) 

where 

l-x ^ 
"i n=l 

is the dilogarithm function. 

The Lagrange multipliers are fixed imposing the mean values 



(E) 



N N da a a 

and 



~hyZ= -\nZ + — In (1 - e- (a+A) ) (56) 



Q N 1- e _(a+A:i 

N c = (n c ) = -— In Z = - In — — . (57) 

oX a 1 — e _A 

Then, the thermodynamic entropy is found to be 
A 

Stherm = In Z + a— + XN C 

= — [2 Li 2 (e~ A ) - 2 Li 2 (e- (Q+A) ) + (a + A) ln(l - e" (a+A) ) - A ln(l - e - A )l . 
a 



(5f 



9 It is useful for the thermodynamic limit to introduce as Lagrange multiplier a = j3N, where j3 is the 
inverse "temperature" . 

10 This holds for gauge group U(N). For SU(N), the product would start with j = 2. 



13 



We would like to verify that the semi-classical fermionic entropy ( 1241) correctly repro- 
duces the thermodynamic entropy (1351) . Taking into account (|T5|) and (15U1) . we need to 
find the mean value (cj), 



N 

(^^^-^.(e^-l)- 1 , (59) 



Thus, the fermion density (|T5|) is^] 

u = 1 - e-^'/JV+A) ; ( 60 ) 

and the entropy (1241) becomes, after making the replacement 



inw H ln(l — u) 

u 



s = -lI iY 



= - [Li 2 (e" A ) - Li 2 (e-( Q+A )) - Li 2 (l - e~ A ) + Li 2 (l - e -( Q+A ))l . (61) 
Eq. ( IBTl) is found to agree with ( 1551) by means of the identity [31] 

Li 2 (x) +Li 2 (l -x) = -ln(l -x)lnx + — . (62) 

6 

5.2 Grand canonical superstar ensemble 

The superstar geometry corresponds to a triangular limit curve with constant slope 
X'(Y) = (cj). From ( 1591) we see that this implies a = 0. For this value of a, the 
expressions for the thermodynamic quantities simplify. Defining 

cu = (e x -l)- 1 , (63) 

the mean energy ( 1561) and number of columns ( 1571) become simply 

A = |iV 2 , N c = ujN , (64) 



and the entropy (16T1) reduces to 



S = N[(l + u) ln(l + w) - oj Inu] 

= V2A [(to 1 ' 2 + lu- 1 / 2 ) ln(l + w) - lu 1/2 In w] 



(65) 



It is interesting to note that, for fixed energy A, the entropy is invariant under lo — > 
exchanging rows with columns. The energy and entropy are, as expected, in agreement 
with the superstar formulae (1341) and 



11 For any smooth function f(x), we have (f(x)) = f((x)) + • • • , where the ellipses denote terms that 
vanish in the thermodynamic limit. 



14 



So far, we have dealt only with the mean quantities, which are independent of the 
choice of ensemble (canonical, grand canonical or micro canonical). Differences between 
the ensembles appear when fluctuations around the mean values are considered, which 
are suppressed in the thermodynamic limit. However, the size of fluctuations is of interest 
for the construction of a stretched horizon around the singularity of the dual geometry. 
This is because the stretched horizon roughly surrounds the space-time region, where 
the underlying microstates significantly differ from each other, and where, therefore, the 
singular solution becomes unreliable. 

Let us, for example, consider the variance of the total number of columns, 

N 

Varn c = (n 2 c ) - (n c ) 2 = ^ (( c * c i) ~ ( c i) ( c j)) ■ ( 66 ) 

It is straightforward to find the correlations 

, ,, ,, , . _(l + w) for i = j, 
(ciCj) - (Cj) {Cj) = { _ (67) 

otherwise. 

Hence, we obtain 

Varn c = Nu{l + u) . (68) 

Let us now estimate the size of a suitable stretched horizon. Setting 6n c = ^Var n c 
and using (jUj), we obtain 

P l _ vMi + -) (69) 



L 2 v/jV * 

This is too large in the large- N limit, because it would imply that the entropy f H2i) asso- 
ciated with the stretched horizon is of order iV 3 / 2 , whereas ( |65|) is of order N. Therefore, 
we conclude that the grand canonical ensemble at infinite temperature contains far too 
many microstates. In the following section, we find a way to improve on this point. 

5.3 Restricted superstar ensemble 

In this section, we introduce an ensemble that exhibits the mean values of the grand canon- 
ical superstar ensemble, while containing smaller fluctuations. We shall call it henceforth 
the restricted superstar ensemble. The main idea is to incorporate the condition on the 
number of columns directly into the sum over configurations instead of imposing it via a 
Lagrange multiplier. Hence, let us define the ensembld^ 2 ! 

{Cj}* = | Ci = 0,l,2,...;^c, <iV c | . (70) 

In words, the restricted superstar ensemble contains all configurations whose total number 
of columns does not exceed iV c . The choice of "does not exceed" as opposed to "is equal 
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For gauge group SU(N), the additional constraint c\ — has to be applied. 
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to" puts N c on the same footing for the columns as iV for the rows and renders the 
ensemble definition symmetric under N «-> iV c o 

The mean energy is imposed in the canonical way. This means that, in order to 
obtain the superstar limit curve, the restricted superstar ensemble should be considered 
at infinite temperature, so we set a = from the start. This is a fortunate circumstance, 
because expectation values can be computed quite easily. Let us start with the partition 
function 

7V C N c -c N N c ~c N c-2 

Z = E* = E E - E !• ( 71 ) 

{cj}* cjv=0cjv_i=0 ci=0 

The sums can be performed by making use of the formula [32] 

Ef'ivr;: 1 )- < 72 > 



k=0 

with the final result 



N + N c 
N 



(73) 



In order to calculate expectation values of products of the c/s, in addition to (I72*|) . 
one also needs the identity 

(« + l)f;W + l)(? + !V (74) 



The results for the one- and two-point correlators are 



(c2) = N C (N + 2N C ) 

\ C J/ (N + l)(N + 2) ' 1 ' 

N C (N C - 1) 
(N + l)(N + 2) 



= TTr 7i\71vr /on for MJ- (77) 



We note that, in contrast to the grand canonical superstar ensemble, the numbers of 
columns of different lengths are now correlated due to the ensemble restriction (ITOl) . 
The mean energy, obtained from ( 1491) and ( 1751) . reads 

A = \NN C . (78) 

This is as expected and agrees with the grand canonical superstar ensemble although, 
now, there was no need to take the large- N limit. To compare the entropy, we consider 
the thermodynamic limit, keeping the slope 

^ = NTT < 79 > 
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Remember that Young diagrams with cjv = are allowed. 
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fixed. We find 

S = InZ « iVm(l + lu) + JV c ln(l + , (80) 

where only the leading term in N has been written. This agrees with (I65I) . as expected. 

We want to verify now whether the estimated size of the stretched horizon determined 
by the magnitude of fluctuations is in agreement with the thermodynamic entropy. There- 
fore, let us consider again the fluctuations of the number of columns. The average number 
of columns is found to be 

^ = WTi = uN ' (81) 

with a variance 

NNJN + N c + 1) N 

where the last equality holds in the large- N limit. Translating this into an estimate for 
the size of the stretched horizon by means of (j47p . we find 



PO 1 /Tr y/u(l+U 

~~ T7v Varra c = — — 



(83) 



L 2 N c N 

This is of the correct order of magnitude in iV and implies an estimated horizon entropy 

ED 



S hov w NtcluVTTuj = itV2K\fuo(l + uo) . (84) 

A slightly smaller estimate can be obtained by realizing that the restricted ensemble 
rOl) does not contain Young diagrams with more than N c columns. That is, in order to 
enclose all mircostate Young diagrams of the restricted ensemble, we can safely place the 
stretched horizon at [c.f. (jTBT) ] 

r 2 H = 2{hN c + hN) = L 4 (uj + ^ + l) , (85) 

whereas the radius of the gray disc is identified with 

rl = 2(h (n c ) + UN) = L A (u + 1) . (86) 



Thus, we obtain 



and 



|| = ^«-2) = £, (87) 
S hoi w Nttuj 3/2 = ncuV2K . (88) 



6 Conclusions and Outlook 

In this paper, we have considered pure and mixed state configurations of the LLM class of 
solutions describing the 1/2-BPS sector of M = 4 super Yang Mills theory. Pure state ge- 
ometries are generated by distributions of concentric rings of droplets, while mixed states 
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are characterized by some density function u(r), which is identified with the semi-classical 
one-fermion phase space density. We have constructed a one-to-one correspondence be- 
tween these density functions and (piecewise) monotonic curves that, in the case of pure 
states, precisely delimit the corresponding Young diagrams. Inspired by the fermionic 
phase space picture, we have proposed a unique entropy formula, eq. (124"1) . and verified 
that it reproduces the thermodynamic entropy of a general grand canonical ensemble of 
Young diagrams. 

Further, we have studied the mean size of fluctuations in thermodynamic ensembles 
of Young diagrams in order to produce an estimate for the size of a stretched horizon 
that encloses the fluctuations of the microstate geometries. Our results show that the 
grand canonical superstar ensemble of pi)] has far too many fluctuations implying that 
the entropy would grow like iV 3 / 2 with the number of fermions, while the thermodynamic 
formula yields a growth like N. To improve upon this, we have introduced a restricted 
superstar ensemble by imposing an upper limit on the number of columns. For this 
ensemble, the estimated entropy scales correctly with iV and differs from the thermody- 
namic entropy by a function of the black hole charge parameter u. As a drawback, we 
observe that our estimates (|84p and (|88|) are not invariant under an exchange of rows 
and columns, keeping the energy fixed. It is conceivable that using the micro canonical 
ensemble (with the restriction n c < N c ) would improve on this point, but that ensemble 
is too difficult to handle explicitly. Furthermore, as the stretched horizon can only yield 
an order-of-magnitude estimate, we would not expect it to reproduce the entropy entirely. 

The superstar belongs to the class of extremal small black holes discussed in the 
introduction. Including higher derivative corrections to SUGRA, it is expected that a 
horizon with finite area is generated. It would be very interesting to obtain a geometric 
entropy, which can be compared with the thermodynamic entropy of the superstar ensem- 
ble. To do this, one could, for example, consider the recently constructed higher derivative 
corrections to five- dimensional SUGRA [33]. Using Sen's entropy function approach in 
conformal supergravity language, similar to [7], seems very promising. Another very in- 
teresting question is whether regular solutions continue to be generated by a distribution 
of droplets, once higher derivative corrections are included. 
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